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Abstract 

We discuss a model for the on-site matrix elements of the sp^d^s* tight-binding hamiltonian of a 
strained diamond or zinc-blende crystal or nanostructure. This model features on-site, off-diagonal 
couplings between the s, p and d orbitals, and is able to reproduce the effects of arbitrary strains on 
the band energies and effective masses in the full Brillouin zone. It introduces only a few additional 
parameters and is free from any ambiguities that might arise from the definition of the macroscopic 
strains as a function of the atomic positions. We apply this model to silicon, germanium and their 
alloys as an illustration. In particular, we make a detailed comparison of tight-binding and ab 
initio data on strained Si, Ge and SiGe. 
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I. INTRODUCTION 



The oncoming limits of conventional downscaling of field-effect transistors have strength- 
ened the need for innovative device architectures.- In this context, the use of mechanical 
strains has become an attractive solution to improve the electrical performances by enhanc- 
ing the carrier mobility.-'^ As a matter of fact, strain engineering techniques such as the 
growth of a contact etch stop layer (CESL)^^ or Si channels strained by SiGe source and 
drain extensions are now widely spread in the semiconductor industry. More generally, the 
electronic properties of strained Sii_a.Ge^ layers grown on Sii^yGey buffers are attracting 
much attention.- These heterostructures, that can be integrated into Si-based electronics 
and photonics, indeed offer the opportunity to tune the band gap of the active layer. 

The modeling of the electrical properties of such devices requires a detailed description of 
the effects of strains on the band structure. Over the past decades, the ab initio methods such 
as the density functional theorySi? (DFT) have provided comprehensive information about 
the deformation potentials of semiconductors.-'^'^'^'^'i^ However, such ab initio methods 
require heavy computational resources and are not, therefore, suitable for the calculation 
of the transport properties of large systems. For that reason, the physics and electronic 
device community is actively developing more efficient semi-empirical approaches, such as 
the k ■ p,—^ the empirical pseudopotentiaU^i^i^ or the tight-bindingi^i^ (TB) methods, 
that can work out the electronic structure of strained semiconductors devices. Among these 
semi-empirical approaches, the TB method has long proved successful in predicting the 
electronic properties of semiconductor nanostructures such as nanocrystals or nanowires. 
The use of an atomic orbitals basis set with interactions limited to a few nearest neighbors 
indeed allows the calculation of the wave functions of million atom systems.— The TB 
method is also well suited to quantum transport calculations,— i^^i^iiS^i^^ and to the atomic 
scale description of, e.g., impurities^^i^S or electron-phonon coupling,^!^ In this respect, the 
first nearest neighbors sp^d^s* model is one of the most accurate and efficient TB description 
of semiconductor materials.— 

The effects of strains are accounted for in TB models through the bond length dependence 
of the nearest neighbor parameters V^y (/i and u being two orbitals on different atoms), which 
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is usually fitted to a power law:— 

v,M = v,Ado) i^-jj , (1) 

where d is the distance between the two atoms in the strained crystal and do is the equi- 
librium distance. Although some hydrostatic and uniaxial deformation potentials can be 
reproduced that way,^^ much better accuracy can be achieved with the introduction of 
strain-dependent on-site parameters.— '^'^'^'^ Indeed, hydrostatic strain shifts the average 
potentiaF^ in the crystal, while uniaxial and shear strains split the p or d orbitals of a 
given atom. In their original sp^d^s* parametrization, Jancu et alM therefore introduced a 
term that lifts the degeneracy between the dyz, dxz, and d^y orbitals under uniaxial (001) 
strain. Jancu and Voisin later generalized this approach to uniaxial (111) strain.— These 
hamiltonians, however, feature the macroscopic strains Eais, whose expression as a function 
of the atomic positions (the basic input of the TB method) is not univocal. Boykin et alr^^ 
therefore introduced position-dependent orbital energies in the sp^d^s* hamiltonian. They 
could reproduce that way the valence band deformation potentials a„ and 6^, but did not 
really improve on d^. This limitation is a consequence of the "diagonal" assumption made in 
that model. Uniaxial (111) strain indeed leaves, for example, the px-, Py and pz orbitals of a 
given atom equivalent. It however couples these orbitals off the diagonal of the hamiltonian. 

In this paper, we discuss a model for the on-site matrix elements of the sp^d^s* TB model, 
based on an explicit expression for the crystal field, assuming that the total potential is the 
sum of central, atomic contributions.— It features off-diagonal couplings between different 
orbitals, and is able to reproduce the effects of arbitrary strains on the band energies and 
effective masses at all relevant k-points. It only involves a few additional parameters, is fully 
consistent with the symmetries of the crystal, and is free from any ambiguity that might 
arise from the introduction of the macroscopic strains Eap in an atomistic description. We 
present this model in section [Tll then discuss its properties in section IIIII Finally, We apply 
this model to silicon, germanium and their alloys, which are the most relevant materials 
for microelectronics, in section IIVI We provide detailed comparisons with ah initio data 
on strained Si, Ge and SiGe, and discuss two important problems: the increase of the 
longitudinal effective mass under shear strains (missing in previous TB models), and the 
description of random alloys. 
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II. MODEL 



In this section, we introduce the model for the on-site matrix elements of the sp^d^s* 
tight-binding hamiltonian. For the sake of simplicity, we focus on a homogeneously strained 
diamond or zinc-blende crystal, the application to arbitrary strains and other crystal struc- 
tures being straightforward. We assume that the total potential in the crystal is the sum of 
central, atomic contributions z^i(|r — Rj|) (sublattice 1) and z^2(|i" — Ril) (sublattice 2), R, 
being the atomic positions. In a first nearest neighbor (NN) approximation, the potential 
experienced by the orbitals of atom i on sublattice 1 is therefore: 



This potential shifts the energy of the orbitals and couples them one to each other in the 
strained crystal. In particular, z/(r) might lift the degeneracy between the p or between 
the d orbitals of the atom. Our model is actually based on a first-order expansion of the 
on-site matrix elements of the potential z/(r) as a function of the atomic positions. In the 
following, we calculate the on-site hamiltonian of the p orbitals of sublattice 1 as an example 
(parapraph III Ap . We then discuss the application to other orbitals and crystal structures 
in paragraph III B[ 

A. Case of p orbitals 

Let p^, pI and pi be the p orbitals of atom i, and: 



NN 




(2) 
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where pi and are the p orbitals ahgned (cr) or orthogonal (vr) to the bond axis Rjj 
Rj — Rj. Slater-Koster relations easily yield:— i^i^ 



NN 



NN 



^ NN 



NN 



(4) 



where = x ■ Ry/dy is the cosine director along x. This expression has been arranged 
so that the last (angular) term of the third line is zero in the unstrained material or under 
hydrostatic pressure (where ~ 4/3 whatever the orientation of the crystal with 

respect to the principal axes). We next expand V2{dij) and V2{dij) in powers of dij — do: 



V^id^j] 
V^{d,,] 



V^{do) + -a: 
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We hence get: 
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do 



do 



where^ ap = (a^ + 2ap/3, = l^a^'l^^o) - ^''(4) and p^p'' = 3{a^ - al)/A. The first 
line of Eq. ([6]) is part of the unstrained p orbital energy E^. The second line is actually 
proportional (to first-order in the c/^'s) to the hydrostatic strain, i.e. proportional to the 
relative variation of the volume of the unit cell (also see paragraph IIIip . We thus define 
for convenience: 

AO on o NN , , 



(1) 



^^0 



do 
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where flo is the unstrained volume of the unit cell. The orbital energy therefore reads 
with these assumptions: 

NN 



An 



where pp{d) = + p^''{d - do) /do- 

The equations are similar for Ey and Ez, with lij replaced by rriij 



(8) 



y • Rij/dij and 



z ■ Hij/dij respectively. Off-diagonal couplings between the p orbitals can be obtained 



in the same way. Slater-Koster relations^^ yield for example: 



NN 



j 
NN 

^ ^ Pp ( dij ) TTlij lij , 



(9) 



which is also zero under hydrostatic pressure. 

The on-site, p block matrix finally reads in the {Px,Py,Pz} basis set: 



I 



NN 



/2 - I Im 
ml m^ — 



In 

mn 



nl nm n^ — h 



(10) 



where the explicit dependence of dij, lij, mij and nij on the atomic sites i and j has been 
dropped for simplicity. The p orbitals feature a oc hydrostatic correction and a oc /3p 
angular term, whose effects will be discussed in more detail in section IIIII 



B. Case of other orbitals 



The on-site hamiltonians of the s (s*) and d orbitals, as well as the off-diagonal coupling 
matrices between the s, p, d and s* orbitals are given in appendix \M Eq- fIlOP as well as 
Eqs. ( IAip - (IA7p of appendix |X] are valid for both sublattices 1 and 2, possibly with different 
parameters in III-V or II- VI materials. They feature hydrostatic (oc a terms) and/or angular 
terms (oc (3 and oc 7 matrices). The oc (3 matrices are all zero in the unstrained crystal and 
under hydrostatic strain. There are, however, non-zero couplings between the d orbitals 
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[Eq. (1A3P ]. between the s and s* orbitals [Eq. (lA4p ]. and between the p and {dyz,dxz,dxy} 
orbitals [Eq. flA7|) ] if the corresponding 7 parameters are not zero. In particular, the d 
orbitals are not degenerate any more in the unstrained crystal if 'yf^ 7^ (see paragraph lA 21 
of the appendix). This is actually consistent with the symmetry of the zinc-blende lattice, 
but is not, usually, accounted for in TB models. As a matter of fact, lifting the degeneracy 
between the d orbitals does not significantly improve the quality of the TB model in diamond 
or zinc-blende crystals, while it is essential in lower symmetry polytypes such as wurtzite 
materials. 

This model has been checked against an ab initio (DFT) description of silicon based on 
atomic-like orbitals (the SIESTA code^'*). With the single-^-polarized basis set used, the 
self-consistent ab initio Hamiltonian is formally equivalent to a non-orthogonal third nearest 
neighbors sp^d^ TB model. The evolution of the on-site ab initio matrix elements under 
strain compares fairly well with our tight-binding approach (despite the latter being first 
nearest neig hbors only). All /5(°)'s and 7(°)'s (except (3^^) are found negative within SIESTA, 
as expected from simple arguments assuming positive, exponentially decaying radial parts 
for the orbitals. The sign of the a's, P^^^'s and 7^^-''s is, however, expected to be quite 
sensitive to the choice of orbitals.— 

The present model can be applied to other crystal structures and inhomogeneous strains. 
In a wurtzite material for example, the oc jSp and oc Pd or 7^ terms will lift the degeneracy 
between the p and between the d orbitals in the unstrained crystal, as is usually enforced a 
priori in the TB descriptions of these materials.— We will now discuss some properties of 
this model, then its application to silicon, germanium and their alloys. 

III. DISCUSSION 

Eq. ffTUl) and Eqs. flAip - (IA7p directly depend on the atomic coordinates through the 
interatomic distances dij and cosine directors kj, rriij, and Uij. These equations are thus 
free of any amibguities that might arise, e.g., from the definition of the strains 
function of the atomic positions, in particular in inhomogeneous environments like alloys. 
They also account for internal strains at the atomistic level, and should therefore be able 
to reproduce electron-optical phonon couplings. Moreover, this model for the on-site tight- 
binding hamiltonian is consistent with the symmetries of the crystal. In particular, the band 
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structure remains invariant under global rotation of the lattice (since these equations fulfill 
Slater-Koster's relationsi^) , a property which is not easily enforced in models depending 



explicitely on the e^/j's or in the model of Refs. |36| and |37|. In practice, the input atomic 
positions can be calculated using, for example, Keating's^ or Stillinger- Weber force fields.— 
We next discuss the effects of biaxial stress on the p orbitals as an illustration of the 
versatility of this model. In a homogeneously strained crystal, the strained atomic positions 
Rj read as a function of the unstrained coordinates R^: 

'R'i = {i + ^ C-{£yz,£xz,£xy), (H) 

where the + (resp. — ) sign holds for sublattice 1 (resp. sublattice 2), is Kleinman's internal 
strain parameter, a is the lattice parameter, I is the identity matrix and e is the matrix of 
the strains Eap- The internal strain parameter ( describes the motion of one sublattice with 
respect to the other under shear strain. We successively consider the cases of biaxial (001) 
and (111) strains. 



A. The case of biaxial (001) strain 



Let us apply a biaxial stress perpendicular to 2; = [001], and let e^x = ^yy = £\\ and 
£zz = £± be the strains in the crystal. Eqs. ( JTOl) and (ITTl) then yield, to first-order in strains: 

H„ = 



-1 











-1 











2 



(12) 



The first line features the hydrostatic strain 6Q/Q0 = e^x + ^yy + ^zz = '^^\\ It accounts 

for the variation of the average potential in the crystal and shifts the three p orbitals equally. 
As expected, the stress also lifts (second line) the degeneracy between the {px^Py} and the 
Pz orbitals. The splitting between {px.Py} and pz, 6Ep = 8l3p^\e± — e\\)/3, is actually 
proportional to the uniaxial component of the strain tensor, but does not depend on 
The degeneracy between {dyz, dxz} and dxy is likewise lifted for the d orbitals. This model 
thus reproduces the effects of the oc (5ooi term in the parametrizations of Jancu et al.,—^ or 
of the diagonal energy shifts in the parametrization of Boykin et al.—^ 
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B. The case of biaxial (111) strain 



Let us now apply a biaxial stress perpendicular to z' = [111]. The strains in the {x' = 
[110], y' = [112], z' = [111]} axis set are thus e^'x' = £y'y' = £\\ and e^'z' = £±- Eqs. f lTUl) and 
(ITT]) then yield, to first-order in strains: 



1 1 

1 1 
1 1 



(13) 



where /3p^ = (3p {1 + 2() + (3p (1 — C)- As expected, biaxial [111] strain leaves the px, Py and 
Pz (diagonal) energies equivalent. It however couples these orbitals off the diagonal of the 
hamiltonian. The eigenvectors of Hp are indeed: i) the p orbital aligned with [111] {Pz'), with 
energy Ep+ap6Q/Qo + 16Pp^{e±—ei\)/27, and ii) the two degenerate p orbitals perpendicular 
to [111] {{Px' , Py'}) , with energies + ap6fl/flo — 8Pp^{e± — e\\)/27. The splitting between 
these orbitals is again proportional to the uniaxial component of the strain tensor. It also 
depends on the internal strain parameter C, (through (3^). The value of C, used as a reference 
to compute the deformation potentials must therefore be provided with the TB parameters. 
Such off-diagonal couplings between the p (or d) orbitals do not exist in the parametriza- 



tion of Ref. 



36 



. As a consequence the degeneracy between the p and between the d orbitals 



is not lifted by biaxial (111) strain, and the value of dy is the same whether the diagonal 



energy corrections are included or not. (3^'' and (3^ also reproduce the effects of the oc 



(1) 
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However, the effective (Sf'^ is assumed 



38 



and oc vrm terms in the parametrization of Ref. 

to be zero for [111] strain (but not for [001] strain), which makes the model of Ref. 
hardly consistent with an explicit description of the crystal field, even beyond first nearest 
neighbors. 



IV. APPLICATION TO SI, GE AND THEIR ALLOYS 

In this section, we discuss the application of the above model for the on-site matrix 
elements of the TB hamiltonian to silicium, germanium and their alloys. We therefore 
attempted to reproduce the band structure of Si, Ge and of the ordered Sio.sGeo.s alloy with 
a first nearest neighbor, two-center orthogonal sp^d^s* TB model. We used experimental 
data when available and ah initio calculations otherwise as a reference for the optimization 
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of the TB parameters. We first review the ab initio calculations and the optimization process 
in paragraph IIV Al then discuss the TB model of Si, Ge and Sio.sGeo.s in paragraph IIV Bl 
and finally the case of arbitrary SiGe alloys in paragraph IIV G[ 

A. First principle calculations and optimization procedure 

A series of first principle calculations was performed with the ABINIT— i^i^ code on 
Si, Ge, and the ordered Sio.sGeo.s alloy, to set up a reference for the optimization of the 
TB parameters. These calculations are based on the local density approximation (LDA) to 
DFT,-i^ using relativistic Hartwigsen-Goedecker-Hutter pseudo-potentials.— The LDA band 
structure was further corrected with Hedin's GW approximation to the self-energy used as 
a post-DFT scheme.— 1^ In general, the GW band energies were found in good agreement 
with the available experimental data.^- The properties of a large set of strained crystals 
have been computed, including hydrostatic as well as biaxial deformations perpendicular to 
[100], [110], and [111].^^ The biaxial strains were chosen large enough (up to = ±5%) to 
span the whole range of lattice mismatches encountered in epitaxial Sii^ajGca- layers grown 
on relaxed Sii_j,Gej, buffers. The atomic positions within the cell were carefully optimized, 
as they strongly affect the band structure.— 

The TB parameters were fitted to the ab initio (or experimental, when available) band 
structures, effective masses and deformation potentials using global optimization methods^ 
refined with local optimizers.— The least-square convergence of the band structures was 
monitored on a dense set of fc-points in the first Brillouin zone. 

The sp^d^s* model of Si and Ge features 4 on-site energies and a parameters, 14 nearest 
neighbor and Harrison (n^u) parameters, and up to 20 (3 and 7 parameters. However, only 9 
of them appeared to have significant impact on the electronic structure of strained Si and Ge 
around the band gap (see Table HTl for a list). In particular, all 7 parameters and most P^^^'s 
were set to zero. This left 45 parameters in the model, that were optimized in following way: 

1. The 4 on-site energies and 14 nearest neighbors parameters were fitted on the band 
structures of relaxed Si and Ge. 

2. The 4 a's and 14 Harrison parameters were fitted on one positive and one negative 
hydrostatic strain.— 
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3. The 7 /?*^°'''s were fitted on one [100] and one [111] biaxial strain that do not change 
the first nearest neighbor bond lengths {e± ~ — Sey and C ^ 1). 

4. The 14 Harrison, 7 /3(°)'s and 2 P^^^'s were further refined on one [100] and two [111] 
biaxial strains [one with ( = 0.557 (Si) or ( = 0.536 (Ge) and one with C = 0]. 

Steps 2 and 3 ensure a reasonnable starting point for step 4. The resulting parametrization 
was also checked against [110] biaxial strains, and its transferability tested on strained Si/Ge 
films and wires. 

The TB model of the Sio.sGeg.s alloy only involves 7 additional first nearest neighbor 
parameters (since the Si/Ge and Ge/Si interactions are different). The on-site energies and 
on-site strain parameters of the Si and Ge atoms were chosen equal to those of bulk Si and 
Ge respectively. 

The TB parameters of Si, Ge and Sio.sGeo.s are listed in Tables HI UTt and lllll The on-site 
strain parameters have the sign expected from simple considerations about the shape of the 
orbitals, except (3p^\ (3sp^ and (3^^] . We point out that the sign of these three parameters is 
extremely robust; Including the missing 7'^°^'s in the on-site corrections will not, in particular, 
change the picture.— The positive sign seems characteristic of first nearest neighbors 

orthogonal models: The model of Ref. indeed splits the p orbitals the same way as ours; 



38l by the choice of an effective (3p^^ = but different effective 



while this is hidden in Ref. 
/3^°^'s for biaxial [001] and [111] strains (see discussion in paragraph IIII Bp . The reasons are 



twofold: First, the orbitals hidden behind orthogonal TB models are much more complex 
than usually assumed when discussing the sign of the interactions. The radial parts must 
indeed have at least one zero to fulfill (near) orthogonality relations with the neighboring 
atoms. Second, some deformation potentials, such as by and H^, are independent on the 
first (but not on the second) nearest neighbor Harisson parameters {b^, for example depends, 
on j3p^^ and jS'^p only). The on-site parameters will therefore likely renormalize the missing 
long-range interactions beyond their "bare" definition given in section III Al and appendix 
lAl The renormalization of long-range interactions into first nearest neighbor and on-site 
terms is underlying every short-range TB model and is a key of their success. We have 
carefully checked our parametrization in bulk (including properties that were not included 
in the optimization, such as the non linearities of the band edges and the behavior of the 
masses under shear strains discussed in the next paragraph), and tested its transferability 
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FIG. 1: (Color online) Band structure of bulk, unstrained silicon in the sp^d^s* TB and GW 
approximations . 



to random SiGe alloys (paragraph IIV Cp and to a variety of test nanostructures such as 
strained Si/Ge films and wires. This model (as the previous ones) actually shows excellent 
transferability of the bulk physics to the nanostructures. 



B. Results in bulk Si, Ge, and SiGe 



The TB and GW band structures of bulk, unstrained Si, Ge and Sio.sGeo.s are compared 
in Figs. [H [21 and [31 They are in very good agreement one with each other, the difference 
between the TB and GW principal band gaps being < 0.01 eV. The lifting of the degen- 
eracies at, e.g., the X point in SiGe are also well reproduced. The TB conduction band 
effective masses and valence band Luttinger parameters of Si and Ge are given in Table 
IIVI They are compared with the GW and experimental data, and with two other sp^d^s* 
parameterizations.— 1^ 

The TB deformation potentials of the conduction and valence band extrema of Si and 
Ge are listed in Table [Vl and compared with the experimental and LDA data. The TB 
model performs well on all relevant deformation potentials. Also shown are the results 
obtained with the sp^d^s* model and parameterization of Refs. [37| and [68|. The present 
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—4.51331 


—4 86118 eV 


V * 


3 06822 


2 92036 pV 
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-0.77711 


—0 23561 eV 


V 

'ppa 


4.08481 


4 36624 eV 
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2.30042 eV 
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-1.56676 


-1.19386 eV 



TABLE I: Tight-binding parameters of relaxed, bulk Si and Ge (first nearest neighbor, two-center 
orthogonal sp^d^s* model). The notations are those of Slater and Koster.— The valence bands of 
Si and Ge have been aligned at = eV; The on-site energies Eg^ Ep, Ed, and Eg* of Ge must, 
therefore, be shifted by Avbo = 0.68 eV to account for the valence band offset between the two 
materials. Ago is the spin-orbit coupling parameter of the p orbitals. 



model reproduces the uniaxial (111) deformation potentials d^j and significantly better, 
as it is able to account for the on-site couplings between the orbitals under shear strains. 



The hydrostatic valence band deformation potential a^, which contro 
band structure on an absolute energy scale, has been fitted to Ref. Il3 



s the position of the 
= 2.38 eV for Si 
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TABLE II: Harrison and on-site strain parameters of Si and Ge. 
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eV 


^P2dia 


-1.61322 


eV 


^pdcr 





75549 


^1^21" 


2 


57110 


eV 


^P2dlTT 


2.43552 


eV 


^pd-K 


1 


67031 



TABLE III: First nearest neighbor two-center tight-binding parameters of SiGe. The on-site ener- 
gies and on-site strain parameters of the Si and Ge atoms are those of bulk Si and Ge respectively. 
The Harrison parameters are the same for Si/Ge and Ge/Si interactions. The relaxed SiGe bond 
length is do = 2.39792 A. 

and = 2.23 eV for Ge). Accordingly, the unstrained valence band offset between Si and Ge 
has been set to Avbo = 0.68 eV, to reproduce the experimental valence band discontinuity 
between Sii-^Ge^. alloys and Si as best as possible (see Table [T] and paragraph IIVC|) . We 
achieve that way strained valence band offsets Avbo = 0.79 eV on Si [001] and Avbo = 0.28 
eV on Ge [001], within the experimental error bars^ and in-between the theoretical LDA 



values of Refs. 



and 



69l . The unstrained valence band offset and hydrostatic deformation 
potential, which are still somewhat controversial,-'^'^'iiii^'i^ can be tuned by shifting all 
on-site energies and a's. 

The typical behavior of the valence and conduction bands of Si and Ge under biaxial [001], 
[110] and [111] strain is plotted as a function of the in-plane deformation e\\ in Figs. [Hand 
O— As a reference, the lattice mismatch of Si grown on Ge is = 4.18 %, while the lattice 
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Si 


Exp. 


GW 


sp^d^s* TB 


Present model 


i 


0.9163'* 


0.925 


0.702^*; 0.891* 


0.900 


i 


0.1905'* 


0.189 


0.227^^; 0.201* 


0.197 


mp 

L 




1.808 


1.378^*; 3.433* 


2.125 






0.124 


O.iei''; 0.174* 


0.151 


7i 


4.26'*; 4.285'* 


4.54 


4.51^*; 4.15* 


4.22 




4.22^4.340'= 








72 


0.38'*; 0.339'* 


0.33 


0.15'*; 0.26* 


0.37 




0.39''; 0.31'= 








73 


1.56'*; 1.446'* 


1.54 


1.55^*; 1.39* 


1.43 




1.44^1.46'= 








Ge 


Exp. 


GW 


sp^d^s* TB 


Present model 




1.588'^; 1.74*= 


1.626 


1.363'*; 1.584* 


1.594 




0.08152'^; 0.079*= 


0.074 


0.083'*; 0.081* 


0.082 






0.881 


0.655''; 0.701* 


0.837 


i 




0.176 


0.223''; 0.201* 


0.178 


uir 






0.038''; 0.039* 


0.038 


7i 


13.0^; 12.8^^; 


13.54 


13.13''; 13.14* 


12.96 












72 


4.4^; 4.05; 


4.32 


4.01'*; 3.68* 


4.11 




4.20** 








73 


5.3^; 5.5^; 


5.77 


5.75'*; 5.63* 


5.59 



5.56'* 



'*Ref. l6l| 

''I. Bdlslev and P. Lawaetz, as presented in Ref. 



m. 



-^Ref. 
'^Ref. M- 
^Ref. 
sRef. ]m 
^Ref. m 
*Ref. 68. 
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TABLE IV: Effective masses and Luttinger parameters of Si and Ge. 
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r X w L r K X 

FIG. 2: (Color online) Band structure of bulk, unstrained germanium in the sp^d^s* TB and GW 
approximations . 

mismatch of Ge grown on Si is £:|| = —4.01%. Under biaxial [001] strain, the six conduction 
band minima of silicon split in two groups,— the A^^y valleys oriented along [100] and [010], 
and the valleys oriented along [001]. The conduction band edges are almost linear with 
strain, A^^y being the lowest energy valleys when e\\ < 0, and when e\\ > 0. Biaxial [110] 
strain also splits the conduction band minima the same way; the A^, minima are, however, 
markedly non-linear, being the lowest energy valleys for both £:|j < and e\\ ^ 2 %. This 
behavior, which is not accounted for by the simplest deformation potential theories, results 
from the shear strains component e^y = {£± — ^||)/2 (see discussion below). As a matter of 
fact, biaxial [111] strain (which does not split the conduction band minimap-^), also exhibits 
the same non-linear trends [Syz = £xz = £xy = {£± ^^||)/3]. It is worthwile to note that these 
non-linearities have not been specifically targeted in the optimization of the TB parameters. 

In Germanium, a biaxial [001] strain likewise splits the A valleys (but not the L ones). 
The Ax^y valleys are the lowest energy bands for compressive strain ^ 2 %, while the F 
valley falls below the L valleys for tensile strain £:|| ^ 2 %. Germanium then becomes a small, 
direct band gap semiconductor, and even a semi-metal (zero gap) when ^ 4 %. Biaxial 
[111] strain splits the L valleys in two groups, ""^^ the three L(lll) valleys (lowest energy for 
compressive strains) and the L[lll] valley (lowest energy for tensile strains). The band gap 
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w 



K X 



FIG. 3: (Color online) Band structure of bulk, unstrained Sio.sGeo.s in the sp^d^s* TB and GW 
approximations . 

again closes when e\\ ^ 3.5 %. The behavior of germanium under biaxial [110] strain is much 
more complex, A^, the two L(lli) valleys, the two L(lll) valleys and F being successively 
the lowest energy band(s) when going from compressive to tensile strains, with a zero band 
gap for e\\ < — 5 % and €\\ > 4 %. 



These results are in good agreement with the ab initio data 
|4]and[5]). Though the sp^d^s parametrization of Refs. l37| and 



black diamonds^ on Figs. 
68 also shows reasonnable 



agreement with ah initio data for [001] biaxial strain, it notably misses the strong non- 
linearity of the conduction band energy in Si [110] or [111]. As stated previously, this 
non-linearity results from the peculiar behavior of the A valleys under shear strains, as 
further evidenced in Fig. [61 Indeed, the conduction band minima move towards the X 
points with increasing compressive or tensile [111] strain, and finally hang to the latter 
when ey > 4 % or e\\ < —3.3 %}^ The position of the conduction band minima along the 
FX RXIS, clS well as the longitudinal effective mass are plotted in Fig. [71 The longitudinal 
effective mass dramatically increases with leyl, and eventually diverges (quartic conduction 
band dispersion) before decreasing again when the conduction band minima reach the edge 
of the first Brillouin zone. Likewise, the A^ valleys shift to the X point under biaxial [110] 
strain, with a divergence of the longitudinal effective mass at £:|j ~ 2.1 % and ~ —2.4 
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Si 




Exp." 


LDA^ 


sp^d^s*"^ 


Present model 


K 




-2.10 ±0.10 


-2.27 


-1.85 


-2.12 






-4.85 ±0.15 


-4.36 


-5.46 


-4.91 


-A _|_ 1;=;A 

' 3' 


— a^j 


1.50 ±0.30 


1.67 


0.97 


1.43 


' 'u 




8.60 ± 0.40 


8.79 


6.88 


8.70 


'—'-7 ~r Q '—'7/ 


- CLv 




-3.14 


-2.61 


-3.20 








13.85 


3.69 


16.19 


Ge 




Exp.''^ 


LDA^ 


sp^d^s*'^ 


Present model 


bv 




-2.86 ±0.15 


-2.90 


-2.48 


-2.74 


dv 




-5.28 ±0.50 


-6 


-3.74 


-5.09 










-9.54 


-9.01 




- a-v 


-2.00 ±0.50 


-2.86 


-2.85 


-3.19 






16.20 ±0.40 


17 


8.09 


15.39 


■jA _|_ I'^A 
^d "T 3""'" 


- ay 




1.43 


3.50 


1.10 


^A 






10 


6.50 


9.02 



'^Cited by Ref. 
^PreseEf 



y 



'^Refs. 
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TABLE V: Deformation potentials of Si and Ge (eV). 



%.i^iI2iiZi The splitting of the transverse masses under uniaxial (110) strains is also well 
reproduced.— These effects, which were not accounted for by previous TB parametrizations, 
are fundamental for the understanding of the transport properties of strained MOSFETS or 
SiGe nanowire heterostructures.-i^i^ 



C. Results in disordered Sii_a:Ge^ alloys 

The TB method offers the opportunity to describe semiconductor alloys as random dis- 
tributions of atoms instead of virtual crystals, thus allowing, e.g., the investigation of alloy 
disorder scattering. The present TB model is particularly well suited to the such random 
alloys as it does not depend on macroscopic strains that would be ill-defined in a disordered 
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Material 


a (N/m) 


f5 (N/m) 


cii (GPa) 


ci2 (GPa) 


C44 (GPa) 


c 


Si 


48.54 


13.83 


165.8 


63.9 


79.3 


0.557 


Ge 


39.14 


11.81 


131.8 


48.3 


64.1 


0.536 


Sio.sGeo.s 


43.80 


12.81 


148.5 


55.9 


71.6 


0.547 



TABLE VI: The valence force field bond-stretching constant a, bond-bending constant /3, elastic 
constants Cij and internal strain parameter of Si, Ge and Sio.sGeo.s. In the disordered SiGe alloys, 
we choose (3 = 13.31 N/m for Si— Si— Ge and [3 = 12.30 N/m for Ge— Ge— Si pairs of bonds. 

environment. We have therefore computed the band gap energy of bulk Sii.ajGe^; alloys 
modeled as random distributions of Si and Ge atoms in large ~ 65000 atoms supercells (in 
order to reduce the statistical noise). The lattice parameters of these supercells and the 
internal coordinates of the atoms were optimized with Keating's valence force field model.— 
The bond bending and bond strecthing constants of the SiGe alloy are given in Table IVIi 

The calculated band gap energy of the alloy is plotted as a function of the Ge mole 
fraction x in Fig. [8|, and compared with luminescence data.~ The lattice parameter of 
the alloy (computed from the valence force field) is also plotted in the inset, and matches 
Dismukes's law^ a(x) = 5.431 + 0.2x + 0.027x^ A (solid line). The present model predicts a 
crossing between the A and L- valley conduction band minima around x = 0.84, in agreement 
with the experimental data. The bowing of the band gap energy for x < 0.84 is, in particular, 
very well reproduced by the tight-binding calculation. We find that the band gap energy 
of the disordered Sio.sGco.s alloy is only ~ 5 meV lower than the band gap energy of the 
ordered alloy. 

As another illustration, the band gap energy of random SiGe alloys biaxially strained on 
Si [001] or [110] is plotted in Fig. [9l Also shown in the inset of Fig. [9]^a) is the valence 
band discontinuity in Sii_a;Gea; alloys on Si [001]. The band gap decreases much faster with 
the Ge mole fraction than in bulk alloys due to the strains. This decrease is again very well 
reproduced by the TB model, showing its ability to describe random alloys under arbitrary 



strains. The conduction and valence band offsets of disordered Sii^ajGe^, al 
[001] buffers are likewise in agreement with the 30 bands k-p model of Ref. 
crystal approximation. 



oy s on Sii__yGej^ 
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FIG. 5: (Color online) (a) [001], (b) [111] and (c) [110] biaxial strain behavior of bulk germanium.— 
The (red) dots, (blue) squares and (magenta) triangles are the L, F and A valleys respectively. 
The black diamonds are the ah initio data.— 
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FIG. 6: (Color online) Lowest two conduction bands (plotted along the FX axis) of bulk silicon 
under [111] biaxial strain. The conduction band minimum is marked with a dot. Its motion with 
strain is plotted as a dotted gray line. The crosses are the ab initio data (also see Fig. [7|). 
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FIG. 7: (Color online) Conduction band, longitudinal effective mass of bulk silicon under [111] 
biaxial strain. Inset: Position of the conduction band minimum along the FX axis. The black 
diamonds are the ab initio data. 
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10 20 30 40 50 60 70 80 90 100 
Ge mole fraction (%) 

FIG. 8: (Golor online) Band gap energy of random, bulk Sii-^Gea; alloys as a function of the 
Ge mole fraction x. Inset: Lattice parameter of the alloy. The solid line is Dismukes's law — 
a{x) = 5.431 + 0.2x + 0.027x2 A. 

V. CONCLUSION 

We have presented a model for the on-site matrix elements of the sp^d^s* TB hamil- 
tonian of a strained diamond or zinc-blende crystal. This model improves over previous 
parametrizations by including the on-site couplings between the s, p and between the d or- 
bitals of the atoms under uniaxial and shear strains. It is able to reproduce the deformation 
potentials and the dependence of the effective masses on strains at all relevant k-points of 
the first Brillouin zone and is fully consistent with the symmetries of the crystal. We have 
succesfully applied this description to Si, Ge and their alloys. This tight-binding model 
should allow predictive modeling of the electronic properties of strained Si/Ge heterostruc- 
tures, and is numerically efficient enough to be included, e.g., in full-band Monte-Carlo^^ or 
Kubo-Greenwood^^ calculations of the transport properties of semiconductor devices. 
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FIG. 9: (Color online) (a) Band gap energy of random Sii-^Ge^: alloys biaxially strained on Si 
[001], as a function of the Ge mole fraction x. The inset is the TB valence band discontinuity AE^ 



(dots) compared to various experimental sources (diamonds) compiled in Ref. 
energy of random Sii-^Gca; alloys biaxially strained on Si [110]. 
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APPENDIX A: COMPLETE ON-SITE COUPLINGS BETWEEN THE s (s*), p, 
AND d ORBITALS 

1. On-site energy of the s (s*) orbitals 

The on-site energy of the s orbitals reads: 

Hs^E',+as^, (Al) 

where is the unstrained s orbital energy and characterizes the dependence of Hg on 
the hydrostatic strain. The same model applies to the s* orbitals (possibly with a different 
a* coefficient). 



2. On-site couplings between the d orbitals 

The on-site d block matrix H.^ reads in the {dy^, d^z, dxy, dx2_y2, d^^'^^^^} basis set: 





' 3 


—Im 


—In 


mn 




NN 


—Im 


2 1 

— ^ 


—mn 


—In 






—In 


—mn 











mn 


—In 







A- 












-n- + i 



NN 



m'^n^ Imn^ Im^n mnu mnv 
Imn^ I'^n^ f'mn Inu Inv 
Im^n Pmn Pm^ Imu Imv 
mnu Inu Imu u^ uv 
mnv Inv Imv uv v^ 
where is the "bare" d orbital energy (see discussion below), u = — m^)/2 and v 
(3n2 - 1)/(2V3). Like pp{d), (5d{d) and -td{d) can be written: 

.(0) , ^(1) d-dp 



(A2) 



(3a{d) = {dl df) - {dj d^) = /3f + 



do 



7,(d) = 3 « dl) + {dl dt) - 4 « dj) = 7f + li 



(1) d — do 



do 



(A3a) 
(A3b) 
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where , /^^ , 7^ and 7^ are additional TB parameters. 

Hrf is the sum of the bare d orbital energies, of a oc hydrostatic correction, and of 
two oc I3di 7d angular matrices. Like Hp, the diagonal of the oc f3d matrix has been shifted 
[by — 1/3 X^ATAT /^rfl*^)] so as to be zero in the unstrained diamond or zinc-blende crystal. 
The five d orbitals are not, however, degenerate any more as soon as 7^'''' 7^ 0. They 
indeed split in two groups: i) the {dyz,dxz,dxy] orbitals with energy + 47^°''/9, and 
ii) the {dx'2-y-2-,d'j,z2_y.2] orbitals with energies E^ (since, e.g., ^^^m^n^ = 4/9).®^ This is 
consistent with the symmetry of the zinc-blende or diamond lattice, though it is usually not 
accounted for in TB models. 



3. Couplings between the s and s* orbitals 

The on-site matrix element coupling the s and s* orbitals reads: 

NN 

Hss*=Y.^sAd), (A4) 

where 7ss*((i) = {si\v2\s*) = 7^°^ + 7is*^^^- is non zero in the unstrained diamond or 
zinc-blende crystal if 7^°* =7^ 0. 



4. Couplings between the s (s*) and p or d orbitals 



The on-site matrices coupling the s and p/d orbitals read: 



NN 



H 



sp 



I m n 



(A5) 



NN 



mn In Im u V 



(A6) 



where (3sp{d) = {s, {u^lp'^) = + f^i^^ and (3sM = Vs {si \u2\ <) 
Both matrices are zero in the unstrained diamond or zinc-blende crystal. 
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5. Couplings between the p and d orbitals 



The on-site matrix coupling the p and d orbitals reads: 



NN 



H 



pd 



Oram / 



n / — m —-^m 



m I 



V3 



n 



NN 



Imn t^n I'^m lu Iv 
m?"n Imn Im^ mu mv 
mrt^ In'^ Imn nu nv 



(A7) 



where Md) = {pj \iy,\ dj) = + /jj^)^ and = V^{p1 Iz/^l - 2 {p- \v^\ dj) 



(0) , 
Tpd + 7, 

.(0) 



(1) d-do 
pd do 



pd do 

. This matrix is non zero in the unstrained diamond or zinc-blende crystal if 



7;.^ ^ 0. 
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